BAR4R-PUB-08/018 
SLAC-PUB-13230 



Observation of e+e — > p + p near -^/s = 10.58 GeV 

B. Aubert, 1 M. Bona, 1 Y. Karyotakis, 1 J. P. Lees, 1 V. Poireau, 1 E. Prcncipe, 1 X. Prudent, 1 V. Tisserand, 1 
J. Garra Tico, 2 E. Grauges, 2 L. Lopez ab , 3 A. Palano ab , 3 M. Pappagallo ab , 3 G. Eigen, 4 B. Stugu, 4 L. Sun, 4 
G. S. Abrams, 5 M. Battaglia, 5 D. N. Brown, 5 R. N. Cahn, 5 R. G. Jacobsen, 5 L. T. Kerth, 5 Yu. G. Kolomensky, 5 
G. Kukartsev, 5 G. Lynch, 5 I. L. Osipenkov, 5 M. T. Ronan, 5 'E K. Tackmann, 5 T. Tanabe, 5 C. M. Hawkes, 6 

N. Soni, 6 A. T. Watson, 6 H. Koch, 7 T. Schroeder, 7 D. Walker, 8 D. J. Asgeirsson, 9 B. G. Fulsom, 9 
C. Hearty, 9 T. S. Mattison, 9 J. A. McKenna, 9 M. Barrett, 10 A. Khan, 10 L. Teodorescu, 10 V. E. Blinov, 11 
A. D. Bukin, 11 A. R. Buzykaev, 11 V. P. Druzhinin, 11 V. B. Golubev, 11 A. P. Onuchin, 11 S. I. Serednyakov, 11 
Yu. I. Skovpen, 11 E. P. Solodov, 11 K. Yu. Todyshev, 11 M. Bondioli, 12 S. Curry, 12 I. Eschrich, 12 D. Kirkby, 12 
A. J. Lankford, 12 P. Lund, 12 M. Mandelkern, 12 E. C. Martin, 12 D. P. Stoker, 12 S. Abachi, 13 C. Buchanan, 13 
J. W. Gary, 14 F. Liu, 14 O. Long, 14 B. C. Shen, 14 'B G. M. Vitug, 14 Z. Yasin, 14 L. Zhang, 14 V. Sharma, 15 
C. Campagnari, 16 T. M. Hong, 16 D. Kovalskyi, 16 M. A. Mazur, 16 J. D. Richman, 16 T. W. Beck, 17 A. M. Eisner, 17 

C. J. Flacco, 17 C. A. Heusch, 17 J. Kroseberg, 17 W. S. Lockman, 17 T. Schalk, 17 B. A. Schumm, 17 A. Seiden, 17 
L. Wang, 17 M. G. Wilson, 17 L. O. Winstrom, 17 C. H. Cheng, 18 D. A. Doll, 18 B. Echenard, 18 F. Fang, 18 

D. G. Hitlin, 18 I. Narsky, 18 T. Piatenko, 18 F. C. Porter, 18 R. Andreassen, 19 G. Mancinelli, 19 B. T. Meadows, 19 
K. Mishra, 19 M. D. Sokoloff, 19 P. C. Bloom, 20 W. T. Ford, 20 A. Gaz, 20 J. F. Hirschauer, 20 A. Kreisel, 20 
M. Nagel, 20 U. Nauenberg, 20 J. G. Smith, 20 K. A. Ulmer, 20 S. R. Wagner, 20 R. Ayad, 21 '0 A. Soffer, 21 -© 
W. H. Toki, 21 R. J. Wilson, 21 D. D. Altenburg, 22 E. Feltresi, 22 A. Hauke, 22 H. Jasper, 22 M. Karbach, 22 J. Merkel, 22 
A. Petzold, 22 B. Spaan, 22 K. Wacker, 22 M. J. Kobel, 23 W. F. Mader, 23 R. Nogowski, 23 K. R. Schubert, 23 
R. Schwierz, 23 J. E. Sundermann, 23 A. Volk, 23 D. Bernard, 24 G. R. Bonneaud, 24 E. Latour, 24 Ch. Thiebaux, 24 
M. Verderi, 24 P. J. Clark, 25 W. Gradl, 25 S. Playfer, 25 J. E. Watson, 25 M. Andreotti afc , 26 D. Bettoni a , 26 C. Bozzi a , 26 
! R. Calabrese ab , 26 A. Cecchi ab , 26 G. Cibinetto ab , 26 P. Franchini ab , 26 E. Luppi ab , 26 M. Negrini ab , 26 A. Petrella ab , 26 
J> ' L. Piemontese", 26 V. Santoro ab , 26 R. Baldini-Ferroli, 27 A. Calcaterra, 27 R. de Sangro, 27 G. Finocchiaro, 27 

S. Pacetti, 27 P. Patteri, 27 I. M. Peruzzi, 27 ! M. Piccolo, 27 M. Rama, 27 A. Zallo, 27 A. Buzzo a , 28 R. Contri ab , 28 
M. Lo Vetere ab , 28 M. M. Macri a , 28 M. R. Monge afc , 28 S. Passaggio a , 28 C. Patrignani ab , 28 E. Robutti a , 28 
A. Santroni Qb , 28 S. Tosi ab , 28 K. S. Chaisanguanthum, 29 M. Morii, 29 J. Marks, 30 S. Schenk, 30 U. Uwer, 30 
\q ! V. Klose, 31 H. M. Lacker, 31 D. J. Bard, 32 P. D. Dauncey, 32 J. A. Nash, 32 W. Panduro Vazquez, 32 M. Tibbetts, 32 
P. K. Behera, 33 X. Chai, 33 M. J. Charles, 33 U. Mallik, 33 J. Cochran, 34 H. B. Crawley, 34 L. Dong, 34 W. T. Meyer, 34 
S. Prell, 34 E. I. Rosenberg, 34 A. E. Rubin, 34 Y. Y. Gao, 35 A. V. Gritsan, 35 Z. J. Guo, 35 C. K. Lae, 35 A. G. Denig, 36 
M. Fritsch, 36 G. Schott, 36 N. Arnaud, 37 J. Bequilleux, 37 A. D'Orazio, 37 M. Davier, 37 J. Firmino da Costa, 37 
G. Grosdidier, 37 A. Hocker, 37 V. Lepeltier, 37 F. Le Diberder, 37 A. M. Lutz, 37 S. Pruvot, 37 P. Roudeau, 37 
M. H. Schune, 37 J. Serrano, 37 V. Sordini, 37 '0 A. Stocchi, 37 G. Wormser, 37 D. J. Lange, 38 D. M. Wright, 38 
I. Bingham, 39 J. P. Burke, 39 C. A. Chavez, 39 J. R. Fry, 39 E. Gabathuler, 39 R. Garnet, 39 D. E. Hutchcroft, 39 

D. J. Payne, 39 C. Touramanis, 39 A. J. Bevan, 40 C. K. Clarke 40 K. A. George, 40 F. Di Lodovico 40 R. Sacco, 40 
M. Sigamani 40 G. Cowan 41 H. U. Flaecher 41 D. A. Hopkins 41 S. Paramesvaran 41 F. Salvatore, 41 A. C. Wren, 41 

D. N. Brown 42 C. L. Davis, 42 K. E. Alwyn 43 D. S. Bailey, 43 R. J. Barlow, 43 Y. M. Chia 43 C. L. Edgar, 43 
G. D. Laffcrty, 43 T. J. West 43 J. I. Yi 43 J. Anderson 44 C. Chen, 44 A. Jawahery 44 D. A. Roberts 44 G. Simi, 44 
J. M. Tuggle, 44 C. Dallapiccola, 45 X. Li, 45 E. Salvati 45 S. Saremi 45 R. Cowan 46 D. Dujmic, 46 P. H. Fisher 46 
K. Koeneke, 46 G. Sciolla, 46 M. Spitznagel, 46 F. Taylor 46 R. K. Yamamoto, 46 M. Zhao, 46 P. M. Patel 47 
S. H. Robertson 47 A. Lazzaro Qb , 48 V. Lombardo Q 48 F. Palombo ab 48 J. M. Bauer 49 L. Cremaldi 49 
V. Eschenburg, 49 R. Godang 49 '0 R. Kroeger 49 D. A. Sanders, 49 D. J. Summers, 49 H. W. Zhao 49 M. Simard, 50 
P. Taras, 50 F. B. Viaud, 50 H. Nicholson, 51 G. De Nardo Qb , 52 L. Lista a , 52 D. Monorchio ab , 52 G. Onorato ab , 52 
C. Sciacca ab , 52 G. Raven, 53 H. L. Snoek, 53 C. P. Jessop, 54 K. J. Knoepfel, 54 J. M. LoSecco, 54 W. F. Wang, 54 
G. Benelli, 55 L. A. Corwin, 55 K. Honscheid, 55 H. Kagan, 55 R. Kass, 55 J. P. Morris, 55 A. M. Rahimi, 55 
J. J. Regensburger, 55 S. J. Sekula, 55 Q. K. Wong, 55 N. L. Blount, 56 J. Brau, 56 R. Frey, 56 O. Igonkina, 56 
J. A. Kolb, 56 M. Lu, 56 R. Rahmat, 56 N. B. Sinev, 56 D. Strom, 56 J. Strube, 56 E. Torrence, 56 G. Castelli ab , 57 
N. Gagliardi Qb , 57 M. Margoni ab , 57 M. Morandin a , 57 M. Posocco a , 57 M. Rotondo", 57 F. Simonetto ab , 57 R. Stroili ab , 57 



2 



C. Voci Qb , 57 P. del Amo Sanchez, 58 E. Ben-Haim, 58 H. Briand, 58 G. Calderini, 58 J. Chauveau, 58 P. David, 58 
L. Del Buono, 58 O. Hamon, 58 Ph. Leruste, 58 J. Ocariz, 58 A. Perez, 58 J. Prendki, 58 L. Gladney, 59 M. Biasini ab , 60 
R. Covarelli ab , 60 E. Manoni ab , 60 C. Angelini ab , 61 G. Batignani ab , 61 S. Bettarini ab , 61 M. Carpinelli ab , 61 '[ll] 
A. Cervelli ab , 61 F. Forti ab , 61 M. A. Giorgi ab , 61 A. Lusiani ac , 61 G. Marchiori ab , 61 M. Morganti ab , 61 N. Neri ab , 61 

E. Paoloni Qb , 61 G. Rizzo Qb , 61 J. J. Walsh a , 61 J. Biesiada, 62 D. Lopes Pegna, 62 C. Lu, 62 J. Olsen, 62 A. J. S. Smith, 62 
A. V. Telnov, 62 F. Anulli a , 63 E. Baracchini ab , 63 G. Cavoto a , 63 D. del Re ab , 63 E. Di Marco ab , 63 R. Faccini ab , 63 

F. Ferrarotto a , 63 F. Ferroni ab , 63 M. Gaspero ab , 63 P. D. Jackson", 63 L. Li Gioi a , 63 M. A. Mazzoni a , 63 S. Morganti a , 63 
G. Piredda a , 63 F. Polci ab , 63 F. Renga ab , 63 C. Voena a , 63 M. Ebert, 64 T. Hartmann, 64 H. Schroder, 64 R. Waldi, 64 

T. Adye, 65 B. Franek, 65 E. O. Olaiya, 65 W. Roethel, 65 F. F. Wilson, 65 S. Emery, 66 M. Escalier, 66 L. Esteve, 66 

A. Gaidot, 66 S. F. Ganzhur, 66 G. Hamel de Monchenault, 66 W. Kozanecki, 66 G. Vasseur, 66 Ch. Yeche, 66 
M. Zito, 66 X. R. Chen, 67 H. Liu, 67 W. Park, 67 M. V. Purohit, 67 R. M. White, 67 J. R. Wilson, 67 M. T. Allen, 68 

D. Aston, 68 R. Bartoldus, 68 P. Bechtle, 68 J. F. Benitez, 68 R. Cenci, 68 J. P. Coleman, 68 M. R. Convery, 68 
J. C. Dingfelder, 68 J. Dorfan, 68 G. P. Dubois-Felsmann, 68 W. Dunwoodie, 68 R. C. Field, 68 A. M. Gabareen, 68 
S. J. Gowdy, 68 M. T. Graham, 68 P. Grenier, 68 C. Hast, 68 W. R. Innes, 68 J. Kaminski, 68 M. H. Kelsey, 68 H. Kim, 68 

P. Kim, 68 M. L. Kocian, 68 D. W. G. S. Leith, 68 S. Li, 68 B. Lindquist, 68 S. Luitz, 68 V. Luth, 68 H. L. Lynch, 68 
D. B. MacFarlane, 68 H. Marsiske, 68 R. Messner, 68 D. R. Muller, 68 H. Neal, 68 S. Nelson, 68 C. P. O'Grady, 68 I. Ofte, 68 

A. Perazzo, 68 M. Perl, 68 B. N. Ratcliff, 68 A. Roodman, 68 A. A. Salnikov, 68 R. H. Schindler, 68 J. Schwiening, 68 
A. Snyder, 68 D. Su, 68 M. K. Sullivan, 68 K. Suzuki, 68 S. K. Swain, 68 J. M. Thompson, 68 J. Va'vra, 68 A. P. Wagner, 68 
M. Weaver, 68 C. A. West, 68 W. J. Wisniewski, 68 M. Wittgen, 68 D. H. Wright, 68 H. W. Wulsin, 68 A. K. Yarritu, 68 

K. Yi, 68 C. C. Young, 68 V. Ziegler, 68 P. R. Burchat, 69 A. J. Edwards, 69 S. A. Majewski, 69 T. S. Miyashita, 69 
B. A. Petersen, 69 L. Wilden, 69 S. Ahmed, 70 M. S. Alam, 70 J. A. Ernst, 70 B. Pan, 70 M. A. Saeed, 70 S. B. Zain, 70 
S. M. Spanier, 71 B. J. Wogsland, 71 R. Eckmann, 72 J. L. Ritchie, 72 A. M. Ruland, 72 C. J. Schilling, 72 
R. F. Schwitters, 72 B. W. Drummond, 73 J. M. Izen, 73 X. C. Lou, 73 F. Bianchi ab , 74 D. Gamba ab , 74 M. Pelliccioni ab , 74 
M. Bomben ab , 75 L. Bosisio ab , 75 C. Cartaro ab , 75 G. Delia Ricca ab , 75 L. Lanceri ab , 75 L. Vitale ab , 75 V. Azzolini, 76 

N. Lopez-March, 76 F. Martinez- Vidal, 76 D. A. Milanes, 76 A. Oyanguren, 76 J. Albert, 77 Sw. Banerjee, 77 
B. Bhuyan, 77 H. H. F. Choi, 77 K. Hamano, 77 R. Kowalewski, 77 M. J. Lewczuk, 77 I. M. Nugent, 77 J. M. Roney, 77 
R. J. Sobie, 77 T. J. Gershon, 78 P. F. Harrison, 78 J. Hie, 78 T. E. Latham, 78 G. B. Mohanty, 78 H. R. Band, 79 
X. Chen, 79 S. Dasu, 79 K. T. Flood, 79 Y. Pan, 79 M. Pierini, 79 R. Prepost, 79 C. O. Vuosalo, 79 and S. L. Wu 79 

(The BaBAR Collaboration) 

1 Laboratoire de Physique des Particules, IN2P3/CNRS et Universite de Savoie, F-7494-1 Annecy-Le-Vieux, France 
2 Universitat de Barcelona, Facultat de Fisica, Departament ECM, E-08028 Barcelona, Spain 
3 INFN Sezione di Bari a ; Dipartmento di Fisica, Universitd di Barf , 1-10126 Bari, Italy 
* University of Bergen, Institute of Physics, N-5007 Bergen, Norway 
5 Lawrence Berkeley National Laboratory and University of California, Berkeley, California 94720, USA 
6 University of Birmingham, Birmingham, B15 2TT, United Kingdom 
7 Ruhr Universitat Bochum, Institut fur Experimentalphysik 1, D-44780 Bochum, Germany 
8 University of Bristol, Bristol BS8 1TL, United Kingdom 
9 University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1 
10 Brunei University, Uxbridge, Middlesex UBS 3PH, United Kingdom 

11 Budker Institute of Nuclear Physics, Novosibirsk 630090, Russia 

12 University of California at Irvine, Irvine, California 92691, USA 
University of California at Los Angeles, Los Angeles, California 90024, USA 

•* University of California at Riverside, Riverside, California 92521, USA 



15 



University of California at San Diego, La Jolla, California 92093, USA 



16 University of California at Santa Barbara, Santa Barbara, California 93106, USA 
17 University of California at Santa Cruz, Institute for Particle Physics, Santa Cruz, California 95064, USA 
18 California Institute of Technology, Pasadena, California 91125, USA 

19 University of Cincinnati, Cincinnati, Ohio 45221, USA 

20 University of Colorado, Boulder, Colorado 80309, USA 

21 Colorado State University, Fort Collins, Colorado 80523, USA 
22 Technische Universitat Dortmund, Fakultdt Physik, D-44^21 Dortmund, Germany 
28 Technische Universitat Dresden, Institut fur Kern- und Teilchenphysik, D-0 1062 Dresden, Germany 
2Jt Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, F-91128 Palaiseau, France 
2j University of Edinburgh, Edinburgh EH9 3JZ, United Kingdom 
26 INFN Sezione di Ferrara a ; Dipartimento di Fisica, Universitd di Ferrara b , 1-44100 Ferrara, Italy 
27 INFN Laboratori Nazionali di Frascati, 1-00044 Frascati, Italy 



3 



28 INFN Sezione di Genova" ; Dipartirnento di Fisica, Universitd di Genova b , 1-16146 Genova, Italy 

29 Harvard University, Cambridge, Massachusetts 02138, USA 
30 Universitat Heidelberg, Physikalisches Institut, Philosophenweg 12, D-69120 Heidelberg, Germany 
31 Humboldt- Universitdt zu Berlin, Institut fur Physik, Newtonstr. 15, D- 12489 Berlin, Germany 
32 Imperial College London, London, SW7 2AZ, United Kingdom 
33 University of Iowa, Iowa City, Iowa 52242, USA 
34 Iowa State University, Ames, Iowa 50011-3160, USA 
35 Johns Hopkins University, Baltimore, Maryland 21218, USA 
36 Universitat Karlsruhe, Institut fur Experimentelle Kernphysik, D-76021 Karlsruhe, Germany 
37 Laboratoire de I'Accelerateur Lineaire, IN2P3/CNRS et Universite Paris-Sud 11, 
Centre Scientifique d'Orsay, B. P. 34, F-91898 Orsay Cedex, France 
38 Lawrence Livermore National Laboratory, Livermore, California 94550, USA 
39 University of Liverpool, Liverpool L69 7ZE, United Kingdom 
40 Queen Mary, University of London, London, El 4NS, United Kingdom 
41 University of London, Royal Holloway and Bedford New College, Egham, Surrey TW20 OEX, United Kingdom 
42 University of Louisville, Louisville, Kentucky 40292, USA 
43 University of Manchester, Manchester M13 9PL, United Kingdom 

44 University of Maryland, College Park, Maryland 20742, USA 
45 University of Massachusetts, Amherst, Massachusetts 01003, USA 
46 Massachusetts Institute of Technology, Laboratory for Nuclear Science, Cambridge, Massachusetts 02139, USA 
47 McGill University, Montreal, Quebec, Canada H3A 2T8 
48 INFN Sezione di Milano" ; Dipartirnento di Fisica, Universita di Milano b , 1-20133 Milano, Italy 
49 University of Mississippi, University, Mississippi 38677, USA 
50 Universite de Montreal, Physique des Particules, Montreal, Quebec, Canada H3C 337 
51 Mount Holyoke College, South Hadley, Massachusetts 01075, USA 
52 INFN Sezione di Napoli a ; Dipartirnento di Scienze Fisiche, 
Universitd di Napoli Federico IF, 1-80126 Napoli, Italy 
53 NIKHEF, National Institute for Nuclear Physics and High Energy Physics, NL-1009 DB Amsterdam, The Netherlands 

54 University of Notre Dame, Notre Dame, Indiana 46556, USA 
55 Ohio State University, Columbus, Ohio 43210, USA 
56 University of Oregon, Eugene, Oregon 97403, USA 
'' 7 INFN Sezione di Padova a ; Dipartirnento di Fisica, Universita di Padova b , 1-35131 Padova, Italy 
58 Laboratoire de Physique Nucleaire et de Hautes Energies, 
IN2P3/CNRS, Universite Pierre et Marie Curie-Paris6, 
Universite Denis Diderot- Paris7, F-75252 Paris, France 
59 University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA 
60 INFN Sezione di Perugia a ; Dipartirnento di Fisica, Universitd di Perugia b , 1-06100 Perugia, Italy 
61 INFN Sezione di Pisa a ; Dipartirnento di Fisica, 
Universita di Pisa ; Scuola Normale Superiore di Pisa c , 1-56127 Pisa, Italy 
62 Princeton University, Princeton, New Jersey 08544, USA 
63 INFN Sezione di Roma a ; Dipartirnento di Fisica, 
Universita di Roma La Sapienza b , 1-00185 Roma, Italy 
64 Universitat Rostock, D-18051 Rostock, Germany 
65 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
66 DSM/Dapnia, CEA/Saclay, F-91191 Gif-sur-Yvette, France 
67 University of South Carolina, Columbia, South Carolina 29208, USA 
68 Stanford Linear Accelerator Center, Stanford, California 94309, USA 
69 Stanford University, Stanford, California 94305-4060, USA 
70 State University of New York, Albany, New York 12,222, USA 
University of Tennessee, Knoxville, Tennessee 37996, USA 
72 University of Texas at Austin, Austin, Texas 78712, USA 
73 University of Texas at Dallas, Richardson, Texas 75083, USA 
74 INFN Sezione di Torino 11 ; Dipartirnento di Fisica Sperimentale, Universita di Torino b , 1-10125 Torino, Italy 
75 INFN Sezione di Trieste 11 ; Dipartirnento di Fisica, Universitd di Trieste b , 1-34127 Trieste, Italy 
76 IFIC, Universitat de Valencia-CSIC, E-46071 Valencia, Spain 
77 University of Victoria, Victoria, British Columbia, Canada V8W 3P6 
78 Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom 
79 University of Wisconsin, Madison, Wisconsin 53706, USA 
(Dated: June 24, 2008) 

We report the first observation of e + e _ — > p + p~ , in a data sample of 379 fb _1 collected with the 
BABAR detector at the PEP-II e + e~ storage ring at center-of-mass energies near =10.58 GeV. 
We measure a cross section of <r(e + e~^ P + p ) =19.5 ± 1.6(stat) ±3.2(syst) fb. Assuming production 
through single-photon annihilation, there are three independent helicity amplitudes. We measure 
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the ratios of their squared moduli to be |F 00 | 2 : |-Fio| 2 : Fn| 2 = 0.51 ± 0.14(stat) ± 0.07(syst) : 
0.10 ± 0.04(stat) ± O.Ol(syst) : 0.04 ± 0.03(stat) ± O.Ol(syst). The |fbo| 2 result is inconsistent with 
the prediction of 1.0 made by QCD models with a significance of 3.1 standard deviations including 
systematic uncertainties. 

PACS numbers: 13.66.Bc, 13.25.-k, 14.40.Ev 



The exclusive production of J/ipVc and other double- 
charmonium vector-pseudoscalar (VP) pairs in e + e~ col- 
lisions around the T(4:S) mass (y/s « 10.58 GeV) is ob- 
served P, Q at rates approximately ten times larger than 
the rates expected from QCD-based models Various 
theoretical efforts have been made to resolve the discrep- 
ancy Measurements of the process e + e~ — > 4>r] H 
provide information on the e + e~ — > VP process in the 
strange quark sector. Study of the vector-vector (VV) 
process e + e~ — ► p + p~ can provide complementary in- 
formation and test perturbative QCD at the amplitude 
level [6| through investigation of the VV angular distri- 
butions. 

The charge-conjugation (C) even final states p°p° and 
<pp° are produced through e + e~ two-virtual-photon an- 
nihilation (TVPA) @, i, ij. For p+p-, C can be ei- 
ther positive or negative. However, due to the particles' 
charges, the e + e~ — * p + p~ process is unlikely to occur 
via TVPA unless there is either significant final quark 
recombination between the products of the two virtual 
photons or final-state interactions (FSI) Assuming 
production through single-photon annihilation or T(AS) 
decay, the VV final state can be described with three in- 
dependent helicity amplitudes. Any discrepancy between 
the amplitudes predicted by perturbative QCD and the 
experimental measurement might indicate contributions 
from mechanisms such as FSI. Such discrepancies could 
help to better understand the importance of FSI effects 
in B 



VV decays [11 



This analysis uses 343 ftV 1 of e + e~ data collected on 
the T(4S) resonance at 10.58 GeV and 36 fb -1 collected 
40 MeV below (off-resonance) with the JBaBAR detec- 
tor at the SLAC PEP-II asymmetric-energy B factory. 
The BABAR detector is described in detail elsewhere [12J . 
Charged-particle momenta and energy loss are measured 
in the tracking system, which consists of a silicon ver- 
tex tracker (SVT) and a helium-isobutane drift chamber 
(DCH). Electrons and photons are detected in a CsI(Tl) 
calorimeter (EMC). Charged pion candidates are identi- 
fied using likelihoods of specific ionization in the SVT and 
DCH, and of Cherenkov angle and photon counts mea- 
sured in an internally reflecting ring-imaging Cherenkov 
detector. Photons are identified by clusters of energy 
deposited in the EMC that have shapes consistent with 
an electromagnetic shower. The clusters are required to 
be isolated, i.e., geometrically unassociated with charged 
tracks. 

To form the p + p~ final state, we select events with ex- 
actly two well-reconstructed oppositely-charged ^ and 



at least two well- reconstructed 7r° candidates. We require 
the tt^ candidates to have at least 12 DCH hits and a 
laboratory polar angle well within the SVT acceptance 
of 0.41 < 8 < 2.54 radians. The laboratory transverse 
momenta of the tt^ candidates are required to be greater 
than 100 MeV/c. The two charged tracks must both be 
identified as pions. We fit the two charged tracks to a 
common vertex, and require the \ 2 probability to exceed 
0.1%. 

The photon candidates used to reconstruct tt° candi- 
dates are required to have a minimum laboratory energy 
of 100 MeV. The invariant masses of the candidate pho- 
ton pairs are required to be within [0.1, 0.16] GeV/c 2 . 
The masses of 7r° candidates are then constrained to the 
world average value [l3| ■ 

The p ± candidates are formed by combining a tt^ can- 
didate with a 7T° candidate. The production angle 8* is 
defined as the angle between the p + meson direction and 
the incident e~ beam in the e + e~ center-of-mass. The 
p ± helicity angles 8± are defined as the angles in the p^ 
rest frame between the direction of the boost from the 
laboratory frame and the direction of the tt ± . We re- 
quire | cos 8* | < 0.8 and | cos 6*± | < 0.85 because there is 
low signal efficiency outside this fiducial region. 

Figure [T] shows the scatter plot of the invariant mass 
TO 7r + 7r o 7r - 7r o versus the absolute momentum difference 
| Ap\ in the laboratory frame between the 7r + 7r°7r~7r° and 
initial e + e~ systems after requiring the tt^tt masses to 
be less than 1.5 GeV/c 2 . The last requirement eliminates 
a twofold ambiguity in forming the p ± candidates. A few 
per cent of the events have more than one p + or p~ can- 
didate because of multiple 7r°'s. All candidates are kept. 

We accept events from within the rectangular area in- 
dicated in Fig. [T] (|m ff + ff o ff - ff o — y/s\ < 0.28 GeV/c 2 and 
\Ap\ < 0.2 GeV/c). There are a total of 612 candidates 
from 571 events in the Y(4S) and off-resonance samples 
combined. Figure [2] shows the scatter plot of the invari- 
ant masses of ir + ir° and 7r~7r° pairs from the accepted 
candidates. The concentration of candidates in the p + p~ 
mass range indicates p + p~ production. 

We use a two-dimensional maximum likelihood fit to 
extract the signal for e + e _ — > p + p~ . Since the final state 
particle masses are far below the e + e~ collision energy, 
we treat the two-body masses as uncorrelated. The sig- 
nal probability density function (PDF) is constructed as 
a product of two identical one-dimensional PDFs for p + 
and p~ . We use a P-wave relativistic Breit-Wigner for- 
mula to construct the PDF for the p ± resonance: 
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FIG. 1: Scatter plot of m^+^o^-^o us. |Ap| between the 
7T + 7r 7r - 7r° and initial e + e~ systems for the on-resonance 
data. 
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FIG. 2: Scatter plot of m^+^o and m^-^o for the accepted 
events in the combined data. 



F(m) OC 



mT{m) 



(mo 2 — m 2 ) 2 + mo 2 T 2 (m) 



(1) 



r(m) = r 



V m / 



1 + g 2 i? 2 
1 + g 2 i? 2 



where m is the observed pion-pair mass, T is the mass- 
dependent p width, and q is the absolute value of the 
pion candidate momentum in the p candidate rest frame. 
The subscript indicates the value at the central mass 
of the p resonance. R is the Blatt-Weisskopf damping 



3J,ll5| 



radius, which we set to 3 ( GeV/c) -1 

A threshold function q 3 /(l + q 3 a) is used to model the 
background in the p^n^n system, where a is a shape 
parameter. We use a linear function to model the residual 
two-dimensional background: 



where a is a floating parameter, and M = 0.89 GeV/c 2 is 
the midpoint of the tt^tt invariant mass range used in 
the fit. 

In the fit to the data, we fix the mass and width of the 
p ± to the world average values [l3| • The parameters var- 
ied in the fit are: a [a(7r + 7r°)=Q!(7r _ 7r C1 )], o, and the num- 
bers of events for the four components: p + p" , p + 7r~7r°, 
p~7r + 7r°, and the residual background. The mass pro- 
jections on m n +^o and m^-^o from the two-dimensional 
fit are shown in Fig. [3l The extracted number of p + p~ 
signal events is 357 ± 29, with 329 ± 25 in the T(45) res- 
onance sample and 31 ± 14 in the off-resonance sample. 

Assuming that the p + p~ is produced through a J PC — 
1 object (a single-photon or T(4S)), and that C and 
parity P are conserved, there are three independent com- 
plex helicity amplitudes, Foo, Fio, and Fn, where the 
indices indicate the helicities of the p mesons. Fio — 
F±w = F 0± i, Fn = F_i_!, and Fi_i = F_ n = due 
to angular momentum conservation [161]. The angular 
distribution of p + p~ decay products can be expressed 
as: 



dN 



d cos 9* d cos 9+ d(p + d cos 9- dip- 



<x\A +l \ 2 + \A-^ 



A±i = sin 9* cos 9+ cos 9- | Foo| 

+ sin 9* sin 9+ sin 6L cos(^+ + ip_ ) \F U \e lVl1 

+- sin 9+ cos 0- (±(1 =F cos 9*)e lip + 

T(l±cosr)e-' i¥, + )|Fio|e lVl ° 
1 



■ cos 9 + sin6L(±(l =p cos6>*)e 
A 

T(l±cosP)e-*-)\F w \e*> 



(3) 



where tp± is the azimuthal angle that corresponds to the 
helicity (polar) angle 9± defined above. In this coordinate 
system, the incoming electron direction has an azimuthal 
angle tp± of zero. The angles (pn and ipio are the strong 
phases of the amplitudes. Due to limited statistics, we 
examine only the projections and thus lose sensitivity to 
these phases. 

The one-dimensional angular distributions are ob- 
tained from Eq. (|3|) by integrating over all other angles. 
When integrating over the full angular range, the results 
are 



dN 



dcos9* 



sin 2 9* | Fool 2 



+/i ( 1 + cos 2 9*)\F w \ 2 + f 2 sin 2 9*\F n 



dN 



occos 2 0±|F OO | 



dcos9± 

+ (/a + / 4 cos 2 0±)|F 1O | 2 + h sin 2 0±|F n | 2 



(4) 



(5) 



Blm^+^o , m^-^a ) = 1 + afm^+^o — M) + a(m 7r - 7r o — M), dN 

(2) d^I 



cx |Foo| 2 + C/6-/7Cos2^ ± )|F 10 | 2 + / 2 |F 11 | 2 , (6) 
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m n + n o (GeV/c 2 ) m^o (GeV/c 2 ) 

FIG. 3: The invariant mass projections a) m n + w a (m n - n a < 1.5 GeV/c 2 ) and b) m n - n o (m n +^o < 1.5 GeV/c 2 ) for accepted 
events in the combined data. The blue-dashed line is the residual linear background, the red dot-dashed line adds a) p + ir~Tv°, 
b) p~-k + tt°, and the blue-dotted line includes both p + -k~tt° and p~-k + -k°. The red solid line adds the signal. 



where the constants /„ are given in the first row of Ta- 
ble D 



TABLE I: Constants in equations [4] \5\ and [6] 



Integrated region 


h 


h 


h 


h 


h 


h 


ft 


Full range 


2 


2 


1 


l 


l 


4 


1 


Fiducial region 


3.15 4.97 


0.77 


1.66 


1.58 


6.44 


3.15 



To determine the amplitude factors, we perform fits 
of Eqs. |U O and [6] to the data. The fits are performed 
in the fiducial region | cos 6** | < 0.8 and |cos#±| < 0.85. 
Limiting the integration to this fiducial region yields the 
constants shown in the second row of the table. 



We use the sPlot [17[ technique to subtract back- 
grounds in the measured angular distributions. This 
technique assigns a weight to each event (sWeight) for 
each category to which it might belong. The sWeights 
are obtained from the 2D fit to the m^-^o versus m n + n o 
distribution. We subdivide cos 9* and cos6± into bins 
and produce an efficiency table from a phase-space-based 
Monte Carlo (MC) simulation. The event weight is given 
by the sWeight divided by the efficiency. 

The background-subtracted and efficiency-corrected 
distributions for cos 6**, tp±, and cos#± are shown in 
Fig. [4j We perform a simultaneous fit of Eqs. 21 [5j 
and [6] to the five angular distributions, assuming there 
are no correlations between the variables. We return to 
the issue of correlations when we discuss systematic un- 
certainties. In the fit, the amplitudes are constrained 
to satisfy |F 00 | 2 + 4|Fi | 2 + 2|Fn| 2 = 1, since there are 
one -Fboj four Fxq and two Fxx amplitude components. 
The free parameters in the fit are |Foo| 2 , l^io| 2 and the 
overall normalization. The value and error of l-Fiil is 



derived from the fit result and its full covariance ma- 
trix using Fu = |(1 - l^ool 2 - 4|Fi | 2 ). The normal- 
ized amplitudes are: |F 00 | 2 : |Fi | 2 : |Fn| 2 = 0.51 ± 
0.14(stat) ± 0.07(syst) : 0.10 ± 0.04(stat) ± O.Ol(syst) : 
0.04 ± 0.03(stat) ± O.Ol(syst). 

To determine the significance of the result and the sys- 
tematic errors in the fitting procedure, we performed fits 
to multiple sets of events generated according to Eq. [3] 
(toy MC). These studies allow us to assess biases that 
arise because of correlations. We find the biases in the fit- 
ted ratios of squared moduli to be less than 0.002, which 
are included in the systematic errors. Most of these bi- 
ases are due to the imperfect MC efficiency corrections 
that result from the coarse bin size of the efficiency table. 
The statistical uncertainties are scaled using the RMS of 
the pull distributions from the toy MC study. The fitter 
underestimates the statistical uncertainties by approxi- 
mately 6%. Other sources of systematic error, such as 
some of those described below for the cross section, have 
little dependence on angle, and thus are expected to be 
relatively small. We neglect them. 

The measured value of |F o| 2 deviates from 1.0, the 
value predicted by QCD models [f|. From the toy MC 
studies, we determine the statistical probability for this 
deviation to be less than 1 in 3000 experiments, corre- 
sponding to 3.4 standard deviations. Including system- 
atic uncertainties, the significance is 3.1 standard devi- 
ations. This suggests that the production may not be 
dominated by single-photon annihilation as naively ex- 
pected. 

The cross section, including radiative corrections, for 
e + e _ — > p + p~ is calculated from 

N 

° Lxex (1 + 5) ' { ' 

where TV is the number of p + p~ signal events extracted 
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from the combined data, L is the integrated luminosity, 
and e is the signal efficiency obtained from MC simu- 
lation that uses the fully differential angular distribu- 
tion derived from the results of the form factor fit. The 
p + p~ signal efficiency in the fiducial region, without ra- 
diative corrections, is estimated to be 15.0%. The cor- 
rection for initial state radiation, 1 + 5, is calculated ac- 
cording to Ref. [IH and has the value 0.775. Assum- 
ing single-photon production, the radiatively corrected 
cross section near ^/s = 10.58 GcV for m p ± < 1.5 
GeV/c 2 and within |cos6>*| < 0.8, |cos0±| < 0.85 is 
8.3 ± 0.7(stat) ± 0.8(syst) fb. Using Eq. [3 we can scale 
the cross section from our acceptance to the full angular 
ranges, which gives 19.5 ± 1.6(stat) ±3.2(syst) fb, where 
the systematic error includes ±1.7 fb due to the effect of 
the uncertainties in the amplitudes on the extrapolation. 

To study the possibility that the observed signal arises 
from T(4S) — ► p + p~ decay, we scale the off-resonance 
signal to the on-resonance luminosity, and subtract it 
from the on-resonance signal. The resulting number 
of events, 35 ± 135, is consistent with zero. The cor- 
responding branching fraction for T(4S*) — » p + p~ is 
(8.1 ± 29.0) x 10~ 7 . The systematic errors, which may 
be estimated from those given below for the cross sec- 
tion, are negligible for this branching fraction measure- 
ment. Restricting possible results to the physical region 
(> 0), the Bayesian 90% confidence level upper limit is 
5.7 x 10~ 6 . 

The systematic uncertainty on the e + e~ — > p + p~ cross 
section, due to uncertainties in the angular distribution 
fit, is estimated by varying the amplitude values. The 
systematic uncertainty from the two-dimensional fit is es- 
timated from the difference in yield obtained by allowing 
the mean and width of the p resonance mass to vary in 
the fit. The systematic uncertainties due to identifi- 
cation, tracking, and tt° efficiency are estimated based on 
measurements from control data samples. The possible 
background from related modes with extra particles is es- 
timated by using extrapolations from four-particle mass 
sidebands. The systematic uncertainties are summarized 
in Table HD 



TABbE II: Systematic uncertainties on the fiducial region 
cross section of e + e _ —> p + p~ ■ 



Source 



Systematic uncertainty % 



Amplitude fit 
Two-dimensional fit 
Particle Identification 
Tracking efficiency 
7T° efficiency 
p + p~ + X feed-down 
buminosity 
Radiative corrections 



5.2 
2.0 
2.3 
1.0 
6.0 
4.9 
2.0 
1.0 




Total 



10.0 



%.8 -0.4 ^ 0.4 0.8 

cose* 

FIG. 4: The background-subtracted (sWeighted) and 
efficiency-corrected a) cos 8+ b) cos 9- c) ip+ d) and e) 
cos#* distributions. The magenta dashed line is the contri- 
bution from the Foo component, the blue dotted line is for 
the Jio component, the blue dot-dashed line is for the Fn 
component, and the red solid line is the total fit function. 



In summary, we have presented the first observation 
of the exclusive production of p + p~ in e + e~ interac- 
tions near y/s = 10.58 GeV and measured the relative 
amplitudes of the three helicity components. Assum- 
ing production through single-photon annihilation, the 
cross section is measured to be a(e + e~~^ p + p~) =19.5 ± 
1.6(stat) ± 3.2(syst) fb. The 90% confidence level up- 
per limit on the branching fraction B(T(AS) — > p + p~) is 
5.7 x 10~ 6 . Our result for the l-Fbol 2 amplitude is incon- 
sistent by 3.1 standard deviations with the predictions of 
QCD models that assume single-photon production, how- 
ever, indicating that other mechanisms such as TVPA 
with FSI may be important. 
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